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ABSTRACT 



We report on gamma-ray observations in the off-pulse window of the Vela pulsar 
PSR B0833— 45, using 1 1 months of survey data from the Fermi Large Area Telescope 
(LAT). This pulsar is located in the 8° diameter Vela supernova remnant, which con- 
tains several regions of non-thermal emission detected in the radio, X-ray and gamma- 
ray bands. The gamma-ray emission detected by the LAT lies within one of these 
regions, the 2° x 3° area south of the pulsar known as Vela-X. The LAT flux is signifi- 
cantly spatially extended with a best-fit radius of 0.88° ±0.12° for an assumed radially 
symmetric uniform disk. The 200 MeV to 20 GeV LAT spectrum of this source is 
well described by a power-law with a spectral index of 2.41 ± 0.09 ±0.15 and integral 
flux above 100 MeV of (4.73 ± 0.63 ± 1.32) x 10~ 7 cm~ 2 s _1 . The first errors repre- 
sent the statistical error on the fit parameters, while the second ones are the systematic 
uncertainties. Detailed morphological and spectral analyses give strong constraints on 
the energetics and magnetic field of the pulsar wind nebula (PWN) system and favor a 
scenario with two distinct electron populations. 

Subject headings: Vela, pulsars, pulsar wind nebula 



1. Introduction 



The Vela pulsar (PSRB0833-45) at a distance of 290 pc (IDodson et all 120031) is one of the 
closest pulsars to Earth and is therefore studied in great detail. Its period of 89 ms and char- 
acteristic age of r c = 11,000 years make it an archetype of the class of adolescent pulsars. As 



with most other pulsars, the Vela pulsar was first detected through radio observations (|Large et al. 
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Kanbach et al. 



(1994) 



1968) and gam ma-ravs ([Thompson! 1 19751). bu t later studied in detail i n the optical ([Wallace et al 
19771). X-rav (IHarnden and Gorensteinl Il973|) and gamma-ray bands (|Kanbach et al.l (|1980i) and 



The pulsar has a spin-down energy loss rate of 7 x 10 36 ergs s 1 with 



the peak electromagnetic power emitted in the GeV gamma- ray band. Indeed, the Vela pulsar is 
the brightest steady astrophysical source for the Fermi-LAT (lAbdo et al.ll2009al) . The gamma-ray 
properties of the pulsar have been studied in detail with the Fermi-LAT, locating the gamma-ray 
emission far out in the magnetosphere close to the last open field-lines. 

Yet ~99% of the pulsar spindown luminosity is not observed as pulsed photon emission and is 
apparently carried away as a magnetized particle wind. Radio and X-ray observations established 
the presence of large scale diffuse e mission surrounding PSRB0833— 45 - thought to be related t o 



the Vela supernova remnant (SNR) ( Dwarakanathlll99ll : lDuncanlll996l : lAschenbach et al 



These radio observations show that the roughly 8° diameter Vela SNR (|Aschenbach et al. 



1995) 



1995) 



co ntains three dist inct central regions of bright diffuse emission, dubbed Vela-X, Vela-Y and Vela- 
Z (|Rishbethll 19581) . The most intense of these, Vela-X, is an extremely bright (~ 1000 Jy) diffuse 
radio structure of size 2° — 3° located close to PSRB0833— 45. Its radio spectral index is sig- 
nificantly harder than those of Vela-Y and Vela-Z, pointing to a young population of non-thermal 
electrons. Indeed, the flat radio spect ral index, the proximity to the Vela pulsar, and the large degree 
of radio polarization in Vela-X led Iweiler & Panagia ( 1980 ) to first suggest that the diffuse radio 
emission is a PWN formed by a relativistic outflow powered by the spin-down of PSR B0833— 45. 
The deceleration of the pulsar-driven wind as it sweeps up ejecta from the supernova explosion 
generates a termination shock at which the particles are pitch-angle scattered and further acceler- 
ated to ultra-relativistic energies. The PWN emission extends across the electr omagnetic spectrum 
in syn chrotron and inverse Compton components from radio to TeV energies (IGaensler and Slane 
2004 PWNe studies can supply information on particle acceleration in shocks, on evolution of 
the pulsar spindown and on the ambient interstellar gas. 

High angular resolution observations of Vela-X in different wavebands showed a rather com- 
plex morphology. X-ray images taken with the Chandra X-ray telescope revealed further de- 
tails (IHelfand et al.ll200l|) : two toroidal arcs of emission, 17" and 30" away from the pulsar, and a 
4' long collimated feature along the pulsar spin axis, which is interpreted as a jet. These structures 
are embedded in an extended nebula located to the south of the Vela pulsar and observed in soft 
X-rays with the ROSAT X-ray telescope. This bright X-ray and radio structure, usually referred 
to as the "cocoon", has an extension of ~ 0.5° x 1.5° and is generally thought to represent PWN 
flow crushed by the passage of the SNR reverse shock. The offset of the cocoon to the south 
of the pulsar is explained by de nse material to the north of PSRB0833— 45 that prevents a sym- 
metric expansion of the PWN (IBlondin et al.ll200l|) . The X-ray spectrum of the cocoon shows a 
ther mal component with a high-energy power-law tail. The detection of VHE gamma-ray emis- 
sion (lAharonian et al.ll2006|) in the cocoon region, albeit at larger angular scales (58' x 43'; yellow 



-6- 



inner contour in Figure Q]), clearly confirmed the notion of a non-thermal particle population in 
this structure. However, these particles do not easily explain the larger and brighter Vela-X radio 
emission in the surrounding "halo" (blue outer contour in Figured)). This led Ide Jager et al.1 (I2008[) 
to suggest a model with two populations of electrons: one at high energies located on the smaller 
cocoon scale, responsible for the X-ray and TeV emission and a second lower energy population 
extending to larger scales and producing the radio flux. These models made a clear prediction that 
the radio-emitting electrons should be visible in the LAT band through inverse Compton scattering 
of the radio-emitting electrons off ambient photon fields. EGRET, the predecesso r of the Fermi 



LAT, was only able to place upper limits on non-pulsed emis s ion fr om this region (IKanbach et al 



1994|) . Recently, using the AGILE satellite, IPellizzoni et al.1 (|2009|) reported the detection of the 
Vela pulsar wind nebula in the energy range from 100 MeV to 3 GeV. 

Here we report on detection of a significant signal in the Vela pulsar off -pulse emission using 
1 1 months of survey observations with the Fermi-LPR. 



2. LAT description and observations 

The LAT is a gamma-ray telescope that detects photons by conversion into electron-positron 
pairs and operates in the energy range between 20 MeV and 300 GeV. It is made of a high- 
resolution converter tracker (direction measurement of the incident gamma-rays), a CsI(Tl) crystal 
calorimeter (ener gy measurement) an d an anti-coincidence detector to identify the background of 
charged particles (lAtwood et al.1 (120091) ). In comparison to EGRET, the LAT has a larger effective 
area (~ 8000 cm 2 on-axis above 1 GeV), a broader field of view (~ 2.4 sr) and a superior angular 
resolution (~ 0.6° 68% containment at 1 GeV for events convertin g in the front section of the 



tracker). Details of the instrumen ts and data process ing are given in lAtwood et al.1 (|2009|) . The 



on-orbit calibration is described in lAbdo et al.l (12009a). 



ip 
l.( 



The following analysis was performed using 1 1 months of data collected starting August 4, 
2008, and extending until July 4, 2009. Only gamma-rays in the Diffuse class events were selected 
(with the tightest background rejection), and from this sample, we excluded those coming from a 
zenith angle larger than 105° to the detector axis because of the possible contamination from Earth 
albedo photons. We have used P6_V3 post-launch instrument response functions (IRFs), that take 
into account pile-up and accidental coincidence effects in the detector subsystems^. 



'See http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/CiceronelATJRFs/IRF_overview.html 
for more details 



Galactic Longitude (deg.) 



Fig. 1. — 61 -GHz WMAP (archival data) radio sky map of Vela-X in galactic coordinates. The 
position of the Vela pulsar is marked with a cross. The blue outer contour shows the region where 
the integral flux densities and the spectral indexes were computed for the radio data. The extraction 
regions for the spectral analysis of the ASCA data are delimited with green boxes. The yellow inner 
line presents the HESS contour at 68% of the peak value. 



3. Timing solution 



The Vela pulsar is the brightest persistent point source in the gamma-ray sky with pulsed 
photons observed up to 25 GeV The study of Vela-X thus requires us to assign a phase to the 
gamma-ray photons and select those in an off-pulse window. Since the Vela pulsar is young and 
exhibits substantial timing irregularities, phase assignment generally requires a contemporary radio 
ephemeris; such a timing model is produced from observations made with the Parkes 64m radio 
telescope. However, Vela is sufficiently bright to be timed directly in the gamma-rays; for this work 
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and that reported in lAbdo et al.1 (|2009d|) . we chose to use a timing model derived directly from LAT 
observations. We used six gamma-ray times of arrival (TOA) covering the commissioning phase 
of the mission (2008 June 25 through August 4) at 5 day intervals and 24 TOAs spaced at 2-week 
intervals during the survey portion of the mission (2008 August 4 through 2009 July 15). The 
TOAs were fitted to a timing model using TEMP02; the RMS residuals o f the TO As with respect 
to the fitted model is 63 /is. More details can be found in lAbdo et al.l (|2009dh . Pulse phases 
were assigned to the LAT data using the F ermi plug-in provid ed by the LAT team and distributed 
with TEMP02. As shown in Figure 1 of lAbdo et al.l (|2009a|) the pulsar emission is quite faint 
in the phase interval <f) = 0.7 — 1.0. We have used this phase interval for both the spectral and 
morphological analysis. 



4. Results 

The spatial and spectral analysis of the gamma-ray e mission was perform ed using two differ 



ent methods. The first is a maximum- likelihood method (|Mattox et al.l l 1996) implemented in the 
Fermi SSC science tools as the "gtlike" code. The second is an analysis tool developed by the LAT 
team called "Sourcelike". In the latter, likelihood fitting is iterated to the data set to simultaneously 
optimise the position and the extension of a source, assuming spatially extended source models and 
taking into account nearby sources as well as Galactic diffuse and isotropic components in the fits. 
Here, we tried both point source and uniform disk models. Sourcelike can also be used to assess 
the Test Statistic (TS) value and to compute the spectra of both extended and point-like sources. In 
this method, the maximum likelihood is performed in independent energy bands, using a region of 
interest whose size is energy dependent: from 15° at 200 MeV to 3.5° at 50 GeV. 

We used the map cube file glljem_v02.fit to model the Galactic diffuse emission together 
with the corresponding tabulated model isotropiciem_v02.txt for the extragalactic diffuse and the 
residual instrument emissions^. Other versions of the Galactic diffuse models, generated by GAL- 
PROP, are also used to assess systematic errors as discussed in section 1431 Nearby so urces in the 
field of view with a statistical significance larger than 5 o are extracted as described in lAbdo et al.1 



(|2009bf) and taken into account in the study. 



2 Available from 



http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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4.1. Morphology 

In the study of the morphology of an extended source, a major requirement is to have the 
best possible angular resolution. Therefore, we decided to restrict our LAT dataset to events with 
energies above 800 MeV, which further reduces the Galactic diffuse background. Figure|2]presents 
the LAT Test Statistic map of off-pulse emission in the Vela region. The Test Statistic (TS) is 
defined as twice the difference between the log-likelihood L\ obtained by fitting a source model 
plus the background model to the data, and the log-likelihood L obtained by fitting the background 
model only, i.e TS = %L\ - L ). This skymap contains the TS value for a point source at each map 
location, thus giving a measure of the statistical significance for the detection of a gamma-ray 
source in excess of the background. Note that the pulsar (cross) is quite faint in this phase interval. 
The skymap shows bright emission south of PSR B0833— 45 with a fainter extension to the east. 
This gamma-ray complex lies within Vela-X; in particular it is contained within the region that 
remains strong at high radio frequencies (denoted by the WMAP flux contours, see discussion). 
An additional source, still unidentified but coincident with the North-Eastern part of the supernova 
remnant Puppis A, is also visible at position (1, b) = (260.30°, —3.16°) with a TS value of 34.6. 
This source is taken into account in the spectral analysis. 

We determined the source extension using Sourcelike with a uniform disk hypothesis (com- 
pared to the point-source hypothesis). The results of the extension fits are summarized in Tabled] 
The difference in TS between the uniform disk and the point-source hypothesis is 47.9 (which 
converts into a significance of ~ 7a for the source extension) for 800 MeV < E < 20 GeV, which 
demonstrates that the source is significantly extended with respect to the LAT point spread function 
(PSF). The fit extension has a radius of 0.88°±0.12°. We support this conclusion in Figure[3l show- 
ing the radial profile for the LAT data above 800 MeV (from the best source location determined 
for a point source fit) and comparing this with the LAT PSF. 

We have also examined the correspondence of the gamma-ray emission with different source 
morphologies by using gtlike with assumed multi-frequency templates. For this exercise we com- 
pared the TS of the point source and uniform model parameters provided by Sourcelike with values 



deriv ed when using morphological templates from the H.E.S.S. gamma-ray excess map (lAharonian et al 



2006h and the WMAP radio images at 61 -GHz (archival data, see section |4~3|) . The resulting Test 



Statistic values obtained from our maximum likelihood fitting are summarized in Tabled Fitting a 
uniform disk to the data using the best location and size provided by Sourcelike improves the TS by 
40.4 in comparison to the point-source hypothesis, comparable to the improvement in TS between 
D and PS models in Table [Q Replacing the disk with spatial template provided by the H.E.S.S. 
observations decreases the TS with respect to the disk hypothesis (ATS = —31.3), implying that 
the LAT emission does not correspond well to the TeV flux. In contrast, using the radio contours as 
spatial template improves the value of the Test Statistic, but only by ATS = +1 1.7. Thus while the 
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best match is with the radio morphology, as expected from the double electron population scenario, 
we cannot (at high significance) rule out a simple disk morphology. 



Model 


Name 


Energy band (GeV) 


K°) 


b(°) 


Radius (°) 


ATS 


Point Source 


PS 


0.8 - 20.0 


263.03 


-3.27 






Disk 


D 


0.8 - 20.0 


263.34 


-3.11 


0.88 ±0.12 


47.9 



Table 1 : Centroid and extension fits to the LAT data for Vela-X using Sourcelike for events with 
energies above 800 MeV. 



Model 


Name 


TS 


Point Source 


PS 


44.0 


Disk 


D 


84.4 


HESS 




53.1 


WMAP61 GHz 




94.0 



Table 2: Comparison of model likelihood fitting results with gtlike for events with energies above 
800 MeV. For each model, we give the name, and the Test Statistic value (TS). 



4.2. Spectral analysis 

The Fermi-LAI spectral points were obtained by dividing the 200 MeV - 20 GeV range into 
7 logarithmically-spaced energy bins and performing a maximum likelihood spectral analysis in 
each interval, assuming a power-law shape for the source. For this analysis we used the uniform 
disk model from Tabled] to represent the gamma-ray emission observed by the LAT, as discussed 
in section l4~Tl Assuming this spatial shape, the gamma-ray source observed by the LAT is detected 
with a significance of 14cr in the 200 MeV - 20 GeV range. The result, renormalized to the total 
phase interval, is presented in Figure @J To determine the integrated gamma-ray flux we fit a 
power-law spectral model to the data in the energy range 200 MeV - 20 GeV with a maximum 
likelihood analysis. This analysis is more reliable than a direct fit to the spectral points of Figure |4] 
since it accounts for Poisson statistics of the data. The spectrum of Vela-X between 200 MeV and 
20 GeV, assuming the uniform disk model from Table [Q is well described by a power-law with 
a spectral index of 2.41 ± 0.09 ± 0.15 and an integral flux above 100 MeV of (4.73 ± 0.63 ± 
1.32) xl0~ 7 cm~ 2 s _1 (renormalized to the full phase interval). This is in agreement with results 
obtained independently using Sourcelike. The first error is statistical, while the second represents 
our estimate of systematic effects as discussed below. No indication of a spectral cut-off at high 
energy can be detected with the current statistics. This result takes into account the gamma-ray 
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Fig. 2.— Test Statistic (TS) map of the PWN Vela-X with side-length 7°, above 800 MeV. Each 
pixel of this image contains the TS value for the assumption of a point-source located at the pixel 
position. WMAP radio contours at 61 GHz (archival data, see section H31) are overlayed as green 
solid lines. The position of the Vela pulsar, PSR B0833— 45, is indicated with a white cross. An 
unidentified source, coincident with the North-Eastern part of Puppis A is visible at position (1, b) 
= (260.30°, -3.16°) 
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Fig. 3. — Radial profile of the LAT data about the best fit position provided by Sourcelike for 
a point source (1, b) = (263.03°, -3.27°) as reported in Table Q] (E>800 MeV). The LAT PSF 
is overlayed as a red solid line for comparison. The background model is presented as a grey 
histogram and the black dots represent the LAT data. 



emission from the source coincident with Puppis A, which was well modeled as a point source 
emitting a power-law of spectral index 1.97 ± 0.16 and integral flux above 100 MeV of (0.43 ± 
0.16)xl0~ 7 cm~ 2 s _1 (statistical errors only). 

As an attempt to estimate the level of pulsed emission in the off-pulse window, we fitted a 
point source at the position of the Vela pulsar, in addition to the uniform disk representing Vela-X. 
We derived an integral flux above 100 MeV of ~ 3 x 10~ 8 cm~ 2 s _1 for the point source, which 
represents ~ 6% of the flux of Vela-X. 

Fitting a point source only at the position of the Vela pulsar, we get a spectrum well described 
by a power-law with a spectral index of 2.98 ± 0.16 and an integral flux above 100 MeV of (1.48 ± 
0.25) x 10~ 7 cm~ 2 s _1 for the phase inter val 0.7 -1.0 (statist ical errors only). This low flux is in 



agreement with the upper limit reported in lAbdo et al.1 (I2009a|) . 



Three different systematic uncertainties can affect the LAT flux estimation. The main system- 
atic at low energy is due to the uncertainty in the Galactic diffuse emission since Vela-X is located 
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only 2° from the Galactic plane in a region of dense molecular clouds. Different versions of the 
Galactic diffuse emission generated by GALPROP were used to estimate this error. The difference 
with the best fit diffuse model is found to be < 6%. By changing the normalization of the Galactic 
diffuse model artificially by ±6%, we estimate this systematic error to be 25% (0.2 - 0.4 GeV), 
14% (0.4 - 0.8 GeV) and <10% (> 0.8 GeV). The second systematic is related to the morphol- 
ogy of the LAT source. The fact that we do not know the true gamma-ray morphology introduces 
another source of error that becomes dominant when the size of the source is larger than the PSF, 
i.e above 600 MeV for the case of Vela-X. Different spatial shapes have been used to estimate this 
systematic error: a disk, a Gaussian and the radio templates. Our estimate of this uncertainty is 
~25% between 600 MeV and 1 GeV and 30% above 1 GeV. The third uncertainty, common to 
every source analyzed with the LAT data, is due to the uncertainties in the effective area. This sys- 
tematic is estimated by using modified instrument response functions (IRFs) whose effective area 
bracket that of our nominal IRE These 'biased' IRFs are defined by envelopes above and below the 
nominal dependence of the effective area with energy by linearly connecting differences of (10%, 
5%, 20%) at log(E) of (2, 2.75, 4) respectively. We combine these various errors in quadrature to 
obtain our best estimate of the total systematic error at each energy and propagate through to the 
fit model parameters. 



4.3. Supporting Multi-wavelength Measurements 



As a means of better understanding the Vela PWN, we compiled and analyzed multi- wavelength 
data corresponding to the longer wavelength synchrotron counterparts of the sub-GeV-peak (halo) 
and TeV-peak (cocoon) Compton emission. Although their morphologies do vary with waveband, 
we have attempted to form the SEDs of the halo and cocoon of Vela-X by using consistent aper- 
tures. This is important in this complex region and has, apparently, not been the practice in some 
previous studies. Vela-X itself has been traditionally studied at low radio frequ encies where the 
spatial resolution is very poor. However, in an 8.4 GHz Parkes image (Fig 2 of lHalesI (|2004|) ). a 
~ 2.5° x 1 .5° region of bright filamentary emission is visible, roughly coincident with the extended 
LAT flux. We examined archival 5-year WMAP sky maps^j and find that this region appears as a 
distinct concentration in the WMAP all-sky images at 23-, 33-, 41-, 61-, and 94-GHz. As the 
resolution increases to higher frequencies it is increasingly separated into eastern and western sub- 
regions, both well south of the Vela pulsar. We measured a flux for each energy band and estimated 
a flux error (dominated by the uncertainty in the background estimation) u sing the r egion defined 
in Figured] This concentration is also clear in the 0.4 GHz all sky maps of lHaslaml (|1982|) . which 



http://lambda.gsfc.nasa.gov 
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Fig. 4. — Spectral energy distribution of Vela-X renormalized to the total phase interval. The LAT 
spectral points are obtained using the maximum likelihood method described in section |4l2l into 7 
logarithmically- spaced energy bins. The statistical errors are shown in blue, while the black lines 
take into account both the statistical and systematic errors as discussed in section I4.2I The red 
dotted line presents the result obtained by fitting a power-law to the data in the 200 MeV-20 GeV 
energy range using a maximum likelihood fit. A 95 % C.L. upper limit is computed when the 
statistical significance is lower than 3 a. 



provide a low-frequency point. The flux measurements are plotted in Figure |5]. We were not able to 
extract a reliable flux estimate from the 8.4 GHz map. We estimate the flux density spectral index 
for this region of Vela-X as a = 0.5 ± 0.05, similar to but steeper than the a = 0.39 ± 0.03 inde x 
measured over 0.03-8GHz for a much larger region covering all of Vela-X (jAlvarez et al.ll200lh . 
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The component measured here is ~ 5 x fainter. Additional mm and IR measurements would be 
very helpful in extending the spectrum and searching for the expected synchrotron peak at ~mm 
wavelengths. 

In the X-ray band many a uthors have estimated the sp ectrum of the cocoon region, start- 
ing with the ROSAT analysis o f IMarkwa rdt & Qgelmanl (|1995|) . More detailed fitting with ASCA 



(|Markwardt & Ogelmanll 19971: iHorns et al.ll2006D showed that the emission must consist of an op- 
tically thin thermal plasma (typically a mekal therma l plasma model) with kT ps 0.3 keV plus a 
power law (T « 2.0 ±0.3) component. More recently. LaMassa et al.1 (|2008|) have analyzed XMM 



data of the bright central portion of the cocoon, and fit a thermal plasma (kT = 0.48Iq:o6 keV ) 
= 2.3 ± 0.04). The 0.2-6.5 keV power-law flux that they find corresponds 
11 ergcm~ 2 s _1 ), when scaled up to the area of the bright H.E.S.S. emission 



plus power law (T 
to 9.5 ± 1.2 x 10~ 
considered here. These authors also fit a hydrogen column density of n# 



1.6 



+0.3 
-0.2 



x 10 



2ll 



cm 



which we shall adopt for our analysis. All of the analyses extending above 2 keV have been forced 
to measure only port ions of the long ~ 1.5° cocoon s tructure. A number of older X-ray (Einstein 



( Har nden et alJll985h: HEAOA-4 ffl 



1996); BeppoSAX dMangano et al 



evin e et al.l ll984)) and soft gamma-ray (OSSE (Ide Jager et al. 



2005 )) observations of the Vela plerion possess a large enough 



field of view to encompass the majority of Vela-X. Yet the spectral extraction regions of these ob- 
servations are centered on the Vela pulsar, such that the bright inner PWN contaminates the low 
surface brightness extended nebula and hardens the net spectrum. We therefore refrain from using 
such archival data as an estimate of the Vela-X spectrum. 

We made a first attempt to improve such measurements by fitting to the combined emission 
in several ASCA GIS2/3 pointings that cover the bulk of the cocoon, as presented in FigureQ] Data 
set 25038000 (76 ks livetime) covered the northern region while data sets 23043000 and 23043010 
(combined livetime 134 ks) cover the southern region. Using XSelect version 2.4, we extract data 
sets from two 20' x 20' regions, one each in the north and south which largely covered this region. 
The large GIS FOV allowed us to select background regions well outside of the cocoon but on 
the same detector. We assumed a fixed absorption % = 1.6 x 10 20 cm -2 and fit a mekal thermal 
plasma plus power law to the combined data. The thermal component is fit with kT = 0.511q q4 
keV (single parameter 90% errors) over all datasets; no significant variation is seen in kT for 
independent fits to the northern and southern regions. To best constrain the power-law component 
we restricted the fit to the 2-10keV range - here separate fits gave T = 1.97too5 m me norm 
and a slightly softer 2.15 ± 0.10 in the south, providing weak evidence for aging of the electron 
population as one moves along the cocoon. Finally for comparison with the H.E.S.S. emission, 
we fit to the combined regions, obtaining an average index of 2.06 ± 0.05 and 2-10 keV flux of 
6.7 ± 0.4 x 10 -11 ergcm~ 2 s _1 (scaling up the flux in our extraction aperture to the area of the 
bright H.E.S.S. emission). This corresponds toa0.3-7keV flux of 1.4 ±0.1 x 10 -10 erg cm -2 s -1 , 
in good agreement with previous estimates. The spectral energy distribution (SED) points from 
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the 2-10 keV power-law portion of this fit are plotted in Figure[5] 

Finally, we wish to check for X-ray emission from the larger halo portion of Vela-X covered 
by the radio/LAT component. This very large region is presently well covered only by the ROSAT 
All-Sky Survey (RASS), which is strongly dominated by the bright thermal emission of the Vela 
SNR, particularly at low energy. To produce a bound on the flux we measured the counts within 
the radio/LAT region in the hard-band 0.5-2.0 keV RASS image, subtracting background from 
appropriate surrounding regions. No significant excess counts were found and we convert the 
upper bound on the flux of a T = 2 power-law component using WebPIMMS, obtaining 2.5 x 
1CT 11 erg cm -2 s _1 . This bound is shown by an arrow in Figured 



5. Discussion 



Different scenarios have been proposed to interpret the multi-wavelength observations of 
Vela-X. Horns et al. (2006) proposed a hadronic model wherein the gamma-ray emission is the re- 
sult of the decay of neutral pions produced in proton-proton collisions in the cocoon. However, this 
model requires a particle density larger than 0.6 cm -3 , which seems disfavored by the recent best 
fit est i mate of therma l partic le density of ~ 0.1 cm -3 using XMM observations (|LaMassa et al 



20081). lLaMassa et al] ([2008) proposed a leptonic model with radio and X-ray emissions resulting 
from synchrotron radiation and gamma-ray emission arising from inverse Compton scattering. In 
this model, the authors need a 3-component broken power-law to describe the electron population 
and adequately fit the data. A model with a single break can also rep roduce the multi-wa velength 
data if a separate electron population produces the radio emission dde Jager et al.l 12008). In this 
case, the morphology of the gamma-ray emission observed by Fermi should be sim ilar to that in 



the ra dio since they are produced by the same electron population. In the model of Ide Jager et al. 



d2008h . the low energy electron component has a total energy of 4 x 10 48 erg, while the X-ray /TeV- 
peak component has a total lepton energy of 2 x 10 46 erg. Both employ a magnetic field of 5 fiG. 

Our new Fermi-LAT spectrum and the improved flux estimates for the radio and X-ray emis- 
sion from the two components of our SED (Figure [5]) allow considerable progress in constraining 
the model parameters. First, the steep LAT spectrum disfavors the hadronic scenario. While the 
VHE gamma-ray data can be adequately fit with gammas from pion decay, neither the ASCA nor 
the LAT data can be accounted for by secondary electrons. We therefore require a three-component 
injection (one hadron and two lepton) in this case, along with a quite high magnetic field in the 
cocoon in order to suppress IC scattering of X -ray emitting electron s from providing the domi- 
nant source of VHE gamma-rays. As noted by Ide Jager et al.1 (|2008), the SED strongly supports 
a two-component leptonic model. We have computed the SEDs from evolving power-law elec- 
tron populations, one each for the X-ray/VHE-peak cocoon and radio/sub-GeV-peak halo. In both 
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regions an exponentially cut-off power law is injected at constant luminosity and evolved for the 
1 1 kyr estimated lifetime of the Vela pulsar, subject to synchrotron and Klein-Nishina adjusted 
Compton losses. We ignore any possible adiabatic losses to the electron population, since these 
are quite uncertain and may, in any case, be offset by the compression from the SNR reverse 
shock. IC seed fields include CMB, far IR ( temperature 25 K, d ensity 0.4 eVcm~ 3 ) and starlight 



(tempera ture 6500 K, den sity 0.4 eV cm 6 (|de Jager et al.1120 08)). reasonable for the the locale of 



Vela-X ( Porter et al.ll2006|) . For each region we vary the magnetic field, power-law cutoff energy, 
power-law index, and total lepton energy; we find the best model fit by minimizing the weighted 
chi-squared statistic between model and data points. For each parameter 90% one-dimensional 
errors are subsequently calculated by varying the best-fit value of the given parameter until chi- 
squared increases by 2.71. The a = 0.5 halo radio spectral index suggests an electron power-law 
index close to the classical p = 2a + 1 = 2. The synchrotron/Compton peak ratio of the cocoon 
implies a B = 4 field, with small uncertainty. In fact we adequately match the SED of both 
components with this field and an E~ 2 spectrum. However for the cocoon region we require a 600 
TeV exponential cut-off and total energy 1.5 x 10 46 erg, while the halo requires a lower 100 GeV 
exponential cutoff and a total energy of 5 x 10 48 erg. The peaks of the cocoon component are 
controlled by the cooling break. The halo population does not cool appreciably during the pulsar 
lifetime and the peak energies are controlled by the exponential cut-off of the injected spectrum. 
The X-ray upper limit on this component is not constraining. Note that we do not require a mid- 
range break in the injected spectrum for either component. 

With so many free parameters, such SED fits are usually illustrative, rather than constrain- 
ing. However, with our new LAT detection and improved low energy measurements we are testing 
the plausible injection spectrum for the Vela-X PWN. We list the parameters determined by chi- 
squared fits to the multi-wavelength data and single-parameter fit errors in Table |3] The cocoon 
emission evidently represents significantly cooled electrons, dominated by relatively recent injec- 
tion of high energy electrons from the pulsar and its termination shock. The halo component, on 
the other hand, represents old electrons - these are easily produced over the lifetime of the pulsar 
for any initial spin period < 60 ms. Although it would be very interesting to push the LAT spec- 
tral measurement to lower energy, where the halo spectrum may peak, this will prove very difficult 
even with more exposure, given the poor low energy PSF On the other hand, extension of the radio 
spectrum through the mm band promises to constrain the high energy cut-off of the halo electron 
spectrum. For the cocoon component, scheduled XMM mapping of this region should provide ap- 
preciable improvement in the spectral measurements of the non-thermal X-rays and may extend to 
low enough energy to probe the synchrotron peak. With such refined constraints we should have a 
quite detailed knowledge of the bulk injection from the pulsar and its termination shock. In turn, 
it may be hoped that this, and similar measurements of other PWNe, will help us understand the 
physics of these relativistic outflows. 
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Component 




E c (eV) 


r 


Etat (erg) 


X 2 /DoF 


Halo 
Cocoon 


qo+0.46 
J - yj -0.38 

1 8n+ - 10 

->- SU -0.08 


l.Oli^ x 10 11 
5.691JS x 10- 


1 Q7+0.02 
L - y/ -0.02 

i goe+o.003 


5.05^ x 10 48 
1.50^;°5 x 10 46 


10.7/9 
57.7/15 



Table 3: Multiwavelength SED fit to the PWN components as seen in Figure [51 
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Fig. 5. — Spectral energy distribution of regions within Vela-X from radio to very high energy 
gamma-rays. Upper panel: Emission from the low energy electron population (halo). WMAP and 
GeV gamma-ray points (this paper) are for the large radio-bright portion of Vela-X. The ROSAT 
upper limit (this paper) on the soft X-ray flux of this region is also shown by an arrow. The Comp- 
ton components from scattering on the CMB (magenta long dashed line), dust emission (magenta 
dashed line) and starlight (magenta dotted line) are shown. Lower panel: Synchrotron and Comp- 
ton emission from the high energy electron population (cocoon). X-ray (ASCA observations, this 
paper) and very high energy gamma-ray (Aharonian et al. 2006) points are also from the cocoon 
region. Only CMB (cyan long dashed line) and dust (cyan dashed line) scattered flux is shown as 
the starlight is Klein-Nishina suppressed. 



